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ABSTRACT 

A quantitative model is presented which describes the ocean surface 
us an ensemble of Hat facets with a normal distribution of slopes. The vari- 
ance of the slope distribution is linearly related to frequency up to 35CJHz and 
constant at higher frequencies. These facets are partially covered with an 
absorbing nonpolarized foam layer. Experimental evidence is presented for 


this model. 
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A MODEL FOR THE MICROWAVE EMIaSIVITY OF 
THE OCEAN’S SURFACE AS A FUNCTION OF WIND SPIH) 

INTRODUCTION 

Hie Scanning Multichannel Microwave Radiometer (SMMR) is a five frequency IU.7. IK, 
21. and 37GHz) dual polarized microwave radiometer which was carried aboard the Nimbu^>-7 and 
Seasat satellites both of which were launched in l<)78. The instrument has an KUcm parabolic dish 
antenna which scans its main beam SO” in azimuth along a conical surface with a 42° cone angle 
and a vertical axis. This provides a constant incidence angle of approximately 50° at the Earth's 
surface for the orbital altitudes of the two spacecTaft (ca 600km Seasat, KOOkm Nimbus). Hie 
spatial resolution is proportional to wavelength and varies from approximately 150km at b.oCillz 
to 25 km at 37GHz. The instrument has been desaibed in detail by GUwrsen and Barath. * Tlie 
purpose of this instrument is to measure sea surface temperature and wind speed at tl.e sea surface 
globally even in the presence of clouds and light rain. 

The SMMR, being a radiometer, measures the upwelling thermal microwave radiation, the in- 
tensity of which is characterized by a brightness temperature. The physical significance of a bright- 
ness temperature is illustrated in Figure I. If microwave radiation with an intensity characterised 
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Figure 1. The Effect of an Absorbing Layer on a Microwave Radiance lixpressed as a 
Brightness Temperature Where Ty Is the Brightness Temperature, 7, fi, and t^ the \b- 
sorption Coefficient, Thickness and Thermodynamic Temperature of the Layer 
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by 1'u is incident on an absorbing (but not scattering or reflecting) layer with an absolute thermo- 
dynamic temperature t^, an absorption coefficient 1 and a thickness 6, the intensity coming out 
the other side is given by: 

Tb • T„e->« + tod . (I) 

riiat is, Tp is given by a term representing the attenuation of the incident radiation and a 
coinplementary radiation term proportional to the absolute temperature of the absorber (Rayleigh- 
Jeans approximation). Note that if T|, to the intensity of the radiation is unchanged; the radia- 
tion is in equilibrium with the absorber. The layer in question could be u section of waveguide, 
an antenna, a radome, or a substantially uniform portion of the atmosphere. For computation 
purpt>ses, the atmosphere is typically treated as many such layers. 

A similar relationship holds for reflection at a surface as is illustrated in Figure 2. If the 
downwelling radiation is given by Tj then the radiation upwelling off the surface Tf is given by 

Tt = RT| + Ft, (2) 

where R is the power reflectivity of the surface, F Is the emissivity of the surface and t, is the 


Tt - RT^ + Et 
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Figure 2. 1 he Effect of Reflection of a Microwave Radiance Expressed as a Brightness Tem- 
perature Off a Surface Such as the Ocean Where Ti Is the Downwelling Brightness Tempera- 
ture, Tt the Upwelling Brightness Temperature and E, R and t,u,f the Emissivity, 

Reflectivity and Thermodynamic Temperature of the Surface Respectively 



uhMilute thtfrinodynamic temperature ut the surface. Consuleration of theniuHlynamic equilihrium 
requires that 

R ♦ L - I (3) 

The abikirption properties of the atmosphere have been discusseii by ( hang and Wilheit.* The 
dominant features in the frequency range ot interest here are a weak water vapor resonance cen- 
tered at 22.23S(iH/ and absorption due to non-raining clouds which is approximately proportional 
to the square of the frequency. The entire problem of radiative transfer in the presence of rain is 
much more complicated as scattering as well as absorption must be considered.-^ The renection 
and emission properties of the ocean surface arc the subject of this paper, l-'igure 3 is a schematic 
superposition of the spectra of all the parameters of interest here. The ordinate is the partial 

3Td 

derivative of the upwelling brightness, temperature with respect to the parameter of interest. — 

oPi 

expressed in arbitrary units. The polarity has been chosen to make the effect positive where it is 
important. The frequencies of the SMMR are marked with the arrows. One can see that these 
are well chosen frequencies for sorting out these effects. There is also some information content 
m the polari/.ation of the brightness temperatures but that is not so easily displayed; it is however 
implicitly exploited in schemes to retrieve the various parameters from the brightness temperature 
measurements.”’**’^’^ Simulations based on these retrieval schemes and measured performance of 
the SMMR instrument"*’^’^ indicate that a measurement accuracy of I.S^C is attainable for the 
sea surface temperature. The lowest frequency, b.bOM/. is used in this retrieval thus the spatial 
resoUilion is limited to roughly 150km. Similarly, the surface wind speed can be extracted from 
the measurements to roughly I m/s accuracy. The lowest frequency used for surface wind speed 
is l0.7CiM/, so approximately ‘)0km spatial resolution is attained. .Atmospheric water can be re- 
trieved with a spatial resolution of bOkm and an accuracy of 0. I^gm/cm^ and 4mg/cm^ for the 
vapor and liquid phases respectively. 
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Figure 3. Schematic Superposition of the Spectra of Various Cicopliysical Parameters. P,. 
The arrows indicate the SMMR frequencies. The signs have been chosen to be positive 
in the fret|uency range of primary importance to the given narameter. 
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In the vjruuis disaisMons"’'*”^’^ of the retrieval !iclieinev. it is obvious that the only sittnificant 
uncertainty tn the rci|uireJ nuulelling is the effect of the wind on the microwave emissivity of the 
oceans surface. The purpose of this paper is to document the model used in the version of the re- 
trieval algorithm^ used at launch of the two spacecraft, to justify this motlel as well as possible and 
to examine its limitations. 

nil MODI L 

It is a straightforward problem to calculate the emissivity of a smooth water surface. 1 he di- 

7 K Q 

electric properties of sea water and saline solutions have Iven discussed by many authors. • • 

We will use values derived from the lane and Saxton^ measurements and expresseil in an analytic 
form by ('liang and W illicit." The formalism for calculating the emissivity for a given view angle 
and polarization is the so called I'resnel relations. The resulting emissivity as a function of view 
angle is shown in l-'igure 4 for a frequency of 10.7011/ and a temperature of 285“K. This smooth 
surface model was also used for calculahng the surface temperature and salinity sensitivity curves 
shown III I'igure 3. I hey were calculalckl for a view angle of 50“ and vertical polarization; the 
results are similar for other angles and for horizontal ptdarization. 

However, the ocean's surface is not a smooth suiface, the wind roughens the surtace, and. if 
It is blowing hard enough, partially covers the surface with foam. It is necessary at this point to 
define more precisely what is meant by surface wind speed. The wind varies with height near the 
surface and tlie details of this variation depend on the temperature difference between the air and 
the sea. ^ ^ When the air is warmer than the ocean (stable) there is more wind shear than for the 
op|x>site (unstable) case for a given wind speed at some reference level. The wind for our purposi* 
is measured near the surface ($ 2CK)m) along with the air and sea temperatures. These data are 
then used to calculate the st>-called friction velocity U* at the sea surface by means of the 
(ardone* ^ model. This is then transformed to 20m height assuming that the sc‘a and air tem- 
perature are equal (neutral stability). This las) step is simply a one-for-oiu transformation to 
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Figure 4. Fmissivity as a Function of View Angle for a Smooth 
Water Surface and for the Ocean Surface with 7 and 14 m/s 
Wind Sneed 
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express the friction velocity in more fumiliar terms. Wherever po«'^ib!e in this pjper the winds 
have been converted to this 20 meter neutral stability wind. 

Cox and Munk^* have quantitatively described the distribution of surface slopes as a function 
of wind speed. They found that the surface slopes were normally distributed with a variance, 
given by 

» 0.003 0048 W (4) 

where ^ is the wind speed in meters/ seconds at 20m height. The factor multiplying W in the 
above equation is slightly different from that in the C'ox and Munk^* paper because the winds 
were measured at I2.5M in their work. These measurements were made at visible wavelengths. 
Much of the roughness they observed is at scales very small compared to microwave wavelengths. 
Tlie iiH)dc! presented here requires only a !r: ction of the Cox and Munk roughness at the longer 
microwave wavelengths. Specifically, the slope variance observed at a given microwave frequency 
is: 


o2(f) = (0.3 + 0.02f(r,H/)o*m f<35CH/ 
o2(f) = ojm f> 35911/ 


(5) 


I'o calculate a rough surface rniissivity from this sIoih.* distribution, one simply averages the 
Fresnel relations^^ over the distribution of surface slopes. In doing so, one implicitly ignores 
surface curvature and all structure comparable to a wavelength and thereby icduces the problem 
to geometric optics. The comparison with obser ations which follow will demonstrate that this 
is a surprisingly good approximation. 

Wind also creates foam on the ocean’s surface. Nordberg, et al. .*^ found a linear increase in 
brightness temperature with wind speed whenever the wind speed exceeded 7 m/s. They were 
viewing directly at the nadir which essentially eliminates the roughness effect leaving foam as the 
most reasonable explanation. In our model, we will treat foam as partially obscuring the surface 
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in j manner imlcpendent of polarization. A non-rcllcclins material partially covering the Mirface 
would have this property as would an absorbing but partially transparent medium with the same 
temperature as the water, bather description alone would be inadequate, but a combination of 
the two descriptions would be closer to reality. The degree to which foam obscures the surface 
is frequency dependent and proportional to the amount by which the wind speed exceeds 7 m/s. 

A reasonable approximation to the available observations of the fraction, K. by which the surfat'C 
relleclivity is reduced by foam is: 


K » a(l - e**^^^“Kw - 7m/s) w>7m/s 
K ■ 0 w <. 7 m/s 


( 6 ) 


where f is the frequency 
and a * 0.(K)6s/m 


f„ » 7.5 GHz 


bmissivities calculated according to this nuxlel for 7 and 1 4 m/s are shown in Figure 4 for 
comparison with the emissivity of a smooth surface. 


SLI’PORriNC. OBSFRVATIONS 


Because the assumed foam model has no polarization character, ucal polarized observations 
of the surface provide a test of the rough surface portion of the model. If one makes the approxi- 
mation that the atmosphere and the surface have the same thermodynamic temperature Tj, then 
it is straightforward to show that for any given view angle 


Ft,(0) 


T,,«^) - T, 

Tv(0) - T, 


Rh(0) 

Rv(0) 


(7) 


Here T„(0) is the horizontal brightness temperature at an angle and K||(0) is the horizontally 
polarized rellectivity. Ty(</) and Ry(0) refer similarly to vertical polarization. Note that because 
Fjj(0) is the ratio of two reflectivities, it is independent of foam cover and thus provides a meas- 
urement of surface roughening. The data from the blectrically Scanned Microwave Radiometer 
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(I SMK) on Nimbua^b (37CiH/. SO” view angle) have been so analyzed and >'ompared with wind 
speeds derived from the operational data buoys. A summary of this comparison is given in 
i'lgure S. The plotted data are for the most part, averages of many observations; a total of 264 
observations are represented. In analyzing the data, it was found that a value of 28S°K for T| 
worked best but that the improvement over any value in the range 280° K was only marginal. 

Uong the model described in the previous section the expected value of Fjg 5 has been calculated: 
the agreement with the observations is striking. .A geometric optics model using the Cox and 
Munk sea surface slope distribution works extremely well at a wavelength of 0.8cm and a view 
angle of 50”. Mollinger*^ has made observations from a fixed platform at frequencies of 1.4. 8.36. 
and l9.34CiHz. He has filtered the data to remove most of the foam effect but application of an 
analysis technique similar to that applied to the Nimbus-6 KSMK data certainly removes the re- 
mainder. These data all can be interpreted in terms of the geometric optic's model but with much 
less slope vananc? than the Cox and .Munk^" values. 

The fiac?ij*is of the Cox and Munk slope variance required to account for the llollinger*^ 
data at 55” incidence angle are plotted in Tigua* 6. The llollinger data are consistent with these 
roughness fractions for all view angles between 0 and 55”. The Nimbus-6 HBMP*"* result is also 
shown in Figure 6. Tliese data form a picture c'cnsistent with tl e roughness required in liquations 
(5) (shown as a solid line). 

The primary available observations relevant to the effect of foam on surfac'e cmissivity are 
from the Bering Sea Expedition (BESEX)*^ and from Cosmos 243.*^ These results, along with a 
plot of dK/3W is given in Figure 7. The observations are difficuit but nevertheless show reasonable 
self-consistency except possibly for the one BESEX point at 37Gllz. 

LIMITATIONS OF THE MODEL 

There are two obvious limitations to this model, the lack of physical optics effects and the 
simplistic treatment of foam. If one calculates the nadir emissivity of the surface according to the 


9 



WIND SPEED 

Figure 5. Observed Value of Fjjj * (Tn - 285)/(Ty - 285) and Wind Speeds. The geometric 
optics curve is based on the Cox and Munk^ ^ distribution of surface slopes. 
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l i^uro (>. IVrcontaui* ol Cox and Munk * “ SU>po Variance Rciiuircd lor a licomotric Optics Model 
lo I xplain the Nimbusf-(> Observations*"* (X) and ilollinger* ^ (•) as a l-'unclion of I'requency 


presi'nt, geonKtric optics, nuKlel there is substantially no change in einissivity through the entire 
0-7 in/s wind s|ved range. Hlume, et al. .**^ have published the obser\ations in |■igua• 8. These 
obserxations at 2.()5Cill/. show a brightness temivratua* increase with increasing wind speed of 
approximately 0.2“K/m/s which corresponds to an apparent einissivity increase of 7 x L'‘^s/m. 
Tven if one allows for the Unite antenna beainwidlh (20..S“) and the increasing contribution of 
rellected sky emission with iiureasing wind induced roughness, only alxnil a quarter of this effect 

I O 

be accounted for. There was not sufficient data in the pa|vr to ivrmit detailed calculations 
of the rellected solar radiation. I'he numbers .ire of the right order of magnitude to account for 
much of the observed effect. In any case, these observation^ sujrgest an up|ier liinii to the physical 
optics contribution of about 5 x IO*-'s/m. 


The treatment of foam here as having neither |M)lari/ation properties nor any view angle 
deivndence is clearly too simple. Williams*** has investigated Mie properties of gelatin stabilized 
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Figure 7. Frequency Dependence of the F'xtent to Which Wind Induced Foam Reduces 
the Microwave Fmissivity of the Ocean Surface 


bubble rafts under laboratory conditions. His results suggest that the increase of emissivity caused 
by foam is caused by the distortion of the meniscus at the foam-water interface which provides a 
gradual transition from the dielectric constant of air to that of water. His a'sults with bubble 
rafts on an aluminum surface suggest that in the bulk of the foam the imaginary part of the index 
of mfrac' on is on the order of 1% that of water. His results form the basis of a calculation of 
foam emissivity to be used as a credibility check on our foam description. Specifically we assume 
a 1 cm thick layer with a complex index of refraction. 

*^foam ~ ^'^water “ 1)/100. (8) 

where is the complex index of refraction of water, to represent the bulk of the foam and a 

linear transition from nfojp, to n^,aier represent the meniscus. The thickness of this linear 
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Figure 8. Brightness Temperature Increase Due to Wind Speed at 
2.65(iH/ from Blumc, et al. 


transition (0.6cm) was adjusted to provide a reasonable fit to the value of f„ useu n Hquation (6). 
The emissivity of this was calculated”^ as a function of frequency and view angle. *'he foam covers 
age efficiency, F|:, was defined in terms of the horizontally polarized emissivity of th»^ foam 
and of smooth water 


Ft 


i.foam I. water 

^tt 'll 


I. water 
‘ It 


(9) 


Foam coverage efficiencies st) calculated are shown in Figure 9 for view angles of 5° and 55°. The 

-f/f 

foam coverage efficiency implied in Fquation (6), ( I - e " ), is shown for comparison. Fxcept for 
some structure which is an artifact of the calculation, (resonances within the I cm foam layer) and 
which certainly would not be observed in natural foam, the agreement is reasonable suggesting 
that ignoring view angle effects in our descriptio.. foam may not be loo serious. 
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FOAM COVERAGE EFFICIENCY 
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Iltc luck of pi>luri/jtion characirr in the foam is more serious. Ky assuming that foam merely 
partiall) obik'urcs the surface we are assuming that the huri/ontal aiul vertical retlectivities are re- 
iluceil m the same pro|Hntion, i.e.: 


i.fvsun 
I V 


(I - I 


foatn 

II 


HI 


i.waiciv 
I V » 



(10) 


In I'lgure 10, the vertical emissivity error made in this assumption is shown for view angles 
of 45“ and 5.5“. IIk* error is less for smalk*r angles and vanishes identically for 0“. I'he resonances 
ohsened in I'igure d aiv ol served in I'igua* 10 and must he similarly discounted. Nevertlhdess 
errors of ahout 0.05 in the vertical emissivity would he expected for lOO'f l\>am cover (wind spee.l 
greater than I70m/s). The ■‘ign of the error is such that the true emissivity in vertical polari/ation 
is somewhat greater than l-()uation (10) would suggest. 


rON( 1 IISION 

A model has been pa'seiiteil for the microwave emissivity of a wind n>ughened. foam covered 
ocean. The roughness |H>rtion of the description is remarkably consistent with obsi’rvations; the 
foam effects show somewhat nu>re scatter. The strength of the foam cover effect at o.b and 10.7 
(ill/ are important parameters in the interpretation of Nimbusr-7 and Seasat SMMK data; the 
strength at higher frei|uencies less so. In comparing the sparx* observations with surface measure- 
ments of temivrature and wind speed, it should be po.ssible to adjust the foam effect at these two 
frei|uencies in order to fine-tune the retrieval algorithm. 
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VERTICAL EMISSIVITY ERROR 



Figure 10. Approximate Frror in Vertical Fmissivity at 100' ? Foam Cover Due 
to Ignoring Polari7.;ition Properties of Foam 
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